Rotating soliton solutions in nonlocal 
nonlinear media 

S. Skupin and M. Grech 

Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, 

Germany 

W. Krolikowski 

OO ' Laser Physics Center, Research School of Physical Sciences and Engineering, Australian 

CD ' National University, Canberra, ACT 0200, Australia 

o 

5_^ . Abstract: We discuss generic properties of rotating nonlinear wave 

O i' solutions, the so called azimuthons, in nonlocal media. Variational methods 

■^ . allow us to derive approximative values for the rotating frequency, which is 

shown to depend crucially on the nonlocal response function. Further on, 
we link families of azimuthons to internal modes of classical non-rotating 
stationary solutions, namely vortex and multipole solitons. This offers an 
exhaustive method to identify azimuthons in a given nonlocal medium. 
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1. Introduction 

Spatial trapping of light in nonlinear media is a result of a balance between diffraction and 
the self-induced nonlinear index change |IT1|2|. As the rate of diffraction is determined by the 
spatial scale of the light beam independently of its dimensionality the stability properties of the 
self-trapped beams critically depend on the particular model of the nonlinear response of the 
medium. Commonly encountered in the context of nonlinear optics the, so called, Kerr nonlin- 
earity with nonlinear index change being proportional to light intensity leads to stable spatial 
solitons only in one transverse dimension. For two dimensional beams the same model does not 
permit stable beam localization. Instead, it predicts either eventual diffraction or catastrophic 
collapse for input power below or above a certain threshold value, respectively ^. Hence, in 
order to ensure the existence of stable self-localized beams in two dimensions different nonlin- 
ear response is necessary. It is known that the saturation of nonlinearity at high intensities is 
sufficient to prevent collapse and provide stabilization of finite beams (see 01 and references 
therein). However, saturable nonlinearity cannot stabilize complex localized structures such as 
vortex beams which are prone to azimuthal instability |5| . It has been shown recently that stabi- 
lization of localized waves is greatly enhanced in the nonlocal nonlinear media. In such media 
nonlinear response in a particular spatial location is typically determined by the wave inten- 
sity in a certain neighborhood of this location. Nonlocality often results from certain transport 
processes such as atomic diffusion i6j, ballistic transport of hot atoms or molecules |7I or heat 
transfer |[8]|9l. It can be also a signature of a long-range interparticle interaction such as in 
nematic liquid crystals ifTOlfTTl [121 . Spatially nonlocal response is also naturally present in 



atomic condensates where it describes a noncontact bosonic interaction lfT3l [T4l [TSJI . Exten- 
sive studies of beam propagation in nonlocal nonlinear media revealed a range of interesting 
features. In particular, it has been shown that nonlocality may affect modulational instability 
of plane waves lfT6l and prevent catastrophic collapse of finite beams ifTTl [TSl as well as sta- 
bilize complex one- and two-dimensional beams including vortices lfT9ll20l 12111221 I23ll24l . 
Recently, it has been shown that nonlocal media can also support stable propagation of the ro- 
tating soli tons, the so called azimuthons 1251 . Azimuthons have been originally proposed in the 
context of local nonUnear media ll26l . They are azimuthally modulated beams with nontrivial 
phase structures exhibiting steady angular rotation upon propagation. They can be considered 
as azimuthally perturbed optical vortices, i.e. beams with singular phase structure 1.271 . Recent 
theoretical studies demonstrated both, stable and unstable evolution of azimuthons Il28ll29l . 
In the former case it has been shown that in a highly nonlocal regimes azimuthons undergo 
structural transformation resulting from the energetic coexistence of solitons of different sym- 
metries. Numerical and analytical studies revealed that the angular velocity of azimuthons is 
governed by two contributions. The linear component, determined solely by the spatial struc- 
ture of the beam (akin to rotation of complex wave structures resulting from beating of their 
constituent modes 1301 ) as well as the nonlinear which is brought about by the nonlinearity l29l . 
In this work we investigate in details rotating nonlocal solitons. We consider various non- 
local models and show how they determine dynamical properties of azimuthons. In particular, 
we will show that both, rotational frequency and intensity profile of the azimuthons can be 
uniquely determined by analyzing Eigenmodes of the linearized version of the corresponding 
nonlocal problem. The paper is organized as follows: In Sec.|2]we introduce the model equa- 
tions, the general ansatz for rotating solitons (azimuthons) and an exact expression for their 
rotation frequency. In Sec. [51 we derive an approximative, variational formula for the rotation 
frequency of the dipole azimuthons. In Sec.lHwe confront our semi-analytical results with ex- 
act numerical solutions of the ensuing nonlinear equations. One of our main findings is that 
azimuthons emerge from internal modes of stationary nonlinear soliton solutions like vortex or 
dipole, which is detailed in Sec.|5l 

2. Rotating Solutions 

We consider physical systems governed by the two-dimensional nonlocal nonlinear 
Schrodinger equation 

d 
/— V/ + Axr+0V' = O. (1) 

oz 

where 9 represents the spatially nonlocal nonlinear response of the medium. Its form depends 

on the details of a particular physical system. 

In the so called Gaussian model of nonlocality, 9 is given as 

d = ^jjc-'^\w{y,z)\^dy, (2) 

where r = xe^ -^ysy denotes the transverse coordinates. 
If 9 is governed by the following diffusion-type equation 

9-^^9 = W\\ (3) 

solving formally Eq. dD in the Fourier space yields 

= ^11 M\r-r'\)\¥ir,z)\'d'^, (4) 

where ^ is the modified Bessel function of the second kind. 



In the following, we will assume that the nonlinear response 9 can be expressed in terms of 
the nonlocal response function R{r) 



9= llR{\r-r'\)\xif{f,zfd^f', 



(5) 



and Eqs. ^ and (|4|i will serve as prominent examples. 

Azimuthons are a straightforward generalization of the usual ansatz for stationary solutions 
(soHtons) 1261 . They represent spatially rotating structures and hence involve an additional 
parameter, the angular frequency Q. 



\lf{r,^,z)^U{r,^-£lzy^ 



(6) 



where U is the complex amplitude function and X the propagation constant. For £2 = 0, az- 
imuthons become ordinary (nonrotating) solitons. The most simple family of azimuthons is the 
one connecting the dipole soliton with the single charged vortex soliton [25J. A single charged 
vortex consists of two equal-amplitude dipole-shaped structures representing real and imagi- 
nary part of U . If these two components differ in amplitudes the resulting structure forms a 
"rotating dipole" azimuthon. If one of the components is zero we deal with the (nonrotating) 
dipole soliton. In the following we will denote the ratio of these two amplitudes by a, which 
also determines the angular modulation depth of the resulting ring-like structure by "1 — a" . 

After inserting the ansatz (|6]l into Eqs. ^ and (|5]), multiplying with U* and d^U* resp., and 
integrating over the transverse coordinates we end up with 



-XL- + QM'+I'+N' = 0. 



(7a) 
(7b) 



This system relates the propagation constant X and the rotation frequency Q. of the azimuthons 
to integrals over their stationary amphtude profiles, namely 
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(7c) 
(7d) 
(7e) 
(7f) 

(7g) 
(7h) 
(7i) 



The first two quantities have straightforward physical meanings, namely "mass" (M) and "an- 
gular momentum" (L-). We can formally solve for the rotation frequency and obtain (for an 
alternative derivation see ||3T1 ) 



D.. 



M{I'+N')-L,{I + N) 
Ll-MM' ■ 



(8) 



Note that this expression is undetermined for a vortex beam (see discussion below). 

While the above formula looks simple its use requires detailed knowledge of the actual so- 
lution which can only be obtained numerically. However, it turns out that it is still possible to 
analyze main features of rotating solitons using simple approximate approach. To this end let 
us consider the "rotating dipole". Obviously, we have Q. = Q for a = (dipole soliton). On the 
other hand, for a = I [vortex soliton y(r)exp(/0 + iX()Z)], we can assume any value for Q. by 
just shifting the propagation constant A = Aq + i2 accordingly (Ao accounts for the propagation 
constant in the non-rotating laboratory frame). However, with respect to the azimuthon in the 
limit a ^ 1, the value of H is fixed. In what follows, we denote this value by n|Q;=i. 

3. Rotation Frequency - Approximate Analysis 

In a first attempt, let us assume we know the radial shape of the vortex soliton V{r). Then, a 
straight forward approximative ansatz for the dipole azimuthons is Il28ll29l 

U{r,(j)- Q.Z) =AV{r) {cos{(j) - Q.z) + /asin(0 - Qz)] , (9) 

where A is an amplitude factor The family of the "rotating dipole" azimuthons has two param- 
eters, e.g. A and Q.. Here, for technical reasons, we prefer to use the mass M and the amplitude 
ratio a. We take the radial shape of the vortex for a given mass Mq and fix the amplitude fac- 
tor A in ansatz (|9|l to A = \/2/{a^ + 1). With this choice the mass of the azimuthon equals 
Mq for all values of a. Then, we can compute all integrals in Eqs. Q as functions of a: 
M = Mo, L- = 2aMo/(a2 + 1), / = k, N ^ [4(1 + a'*)A^cc + ^a^Ncs]/ia^ + 1)^ M' = Mq, 
/' = 2aIo/{a^ + 1), and A^' = 4a{Ncc +A^cs)/(«^ + !)■ The index "0" indicates the value of the 
respective integral for the vortex (i.e., for a = 1), and 

A^cc= nil R{\r-r'\)V^{r')cos^{^')V^{r)cos^{^)cf-r'cf-r (10a) 

A^cs= nil R{\r-r'\)V^{r')cQS^{^')V^{r)s.m^{^)d-r'cfr. (10b) 

Finally, we obtain the following simple expression for the rotation frequency from Eq. ^, 

a^J^m^^M. (u) 

a^ + l Mo 

Let us have a closer look at this approximation. First of all, we have A^cc > A^cs if we as- 
sume R{r) to be monotonically decreasing in r > 0. This immediately implies that the sense 
of the amplitude rotation is opposite to that of the phase rotation. In general, we need to know 
the actual shape of the vortex soliton to compute Eq. ( fTTl i. One possibility would be to use 
exact numerical solutions of the vortex shape V, a time consuming task because V depends 
in nontrivial way on Aq. Moreover, Eq. (fTTI ) is only an approximation to Q. and should not 
require the exact shapes. Luckily, it is possible to compute quite good approximative vor- 
tex solitons using the so called variational approach ll32l . We employ the following ansatz 
V — Arexp(— r^/2(7^), U = V exp(/0), Xj/ = L'^exp(;Aoz) and look for critical points of the La- 
grangian Jff = —XoM + I+N/2 with respect to the variables A and a. In the following, we give 
results for our two exemplary nonlocal nonlinearities. 

For the Gaussian response, Eq. (O, all integrals appearing in the variational approach can be 
evaluated analytically. Hence, we find the width a of our approximate solution as 

AoC7*^ + (-2 + 2Ao) a" + (-8 +4Ao) a'* - Sa^ - 8 = 0, (12) 




1000 



< 500 




500 



20 
15 
10' 



(c) 








^ 




^ 


X 


/- 


^1 


0.5 

0.4 
C 


- 


- - 














200 

^0 


400 




— — I — - 


- — — .- 


_ _ _ 


— -r 


_ _ _ 


_, — _ — _ 



100 200 300 400 500 

Fig. 1. Nonlocal model and variational approach: Solid lines show (a) width CT, (b) ampli- 
tude A, and (c) frequency 0.\a=\ for the Bessel response as functions of the propagation 
constant Ay. Dashed lines represent the same quantities for the Gaussian nonlocal response. 
The inset in (c) shows ii|a=i for the latter response on more appropriate scales. While both, 
width CT and amplitude A, exhibit similar behavior for the two nonlocal models the rotation 
frequency n|(;(=i shows a completely different dependence on Ay and drastic difference in 
magnitude for Gaussian and Bessel responses. 
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and the rotation frequency Q. 
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We show these three quantities as functions of Ao in Fig. [T] (dashed lines). If Aq is large, (7^1 
(nonlocal limit) and we can neglect the first two terms in Eq. (fT4l) . Then we get (7ni = (2/Ao) '' ^, 
^^ni = ■^(<7ni + 1)^'^ and, subsequently, iini = a/ {a? + 1). The immediate conclusion drawn 
from this formula is that dipole azimuthons in the Gaussian nonlocal model do not rotate faster 
than ii ~ 0.5. Moreover, in the nonlocal regime the rotating frequency does not change much 
with the propagation constant A or the mass M of the azimuthon, but is mainly determined by 
the structural parameter a. 

For the Bessel nonlocal response, Eq. (HJi, we have to compute the relevant integrals numeri- 
cally. The results for both, Bessel and Gaussian models are presented in Fig.[T|using solid and 
dashed lines, respectively. It is clear that amplitude A as well as soliton width a, follow similar 
dependencies as a function of the propagation constant Aq [see Fig. [n^-b)]- Totally different 
dependence shows the rotation frequency £1 [see Fig. [He)]. In particular, not only azimuthons 
in the Bessel nonlocal model rotate much faster than those predicted by the Gaussian model but 
also the former exhibit strong sensitivity of £1 on the propagation constant A and mass M. 

To shed light on the origin of the difference in the rotating frequencies of the two models let 
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Fig. 2. (a) Mass M and (b) quantity N^^ — A'cs for Bessel (solid lines) and Gaussian 
(dashed lines) nonlocal responses as functions of the propagation constant Aq. It is ob- 
vious that the difference in the rotating frequency fl|a=i shown in Fig. [TJc) is due to 
A'cc — A'cs [see Eq. dllb l. (c) Spatial Fourier Transform (FT) of the Bessel (solid lines) and 
Gaussian (dashed lines) nonlocal response R as a function of the transverse wave-vector 

kj^ = \ k\ + k^,. The degree of spatial averaging is obviously larger for the Gaussian than 
for the Lorentzian. 



us have a look at the quantities M and A^cc — Nca, which directly determine H | „=! in the present 
approximation [see Eq. (fTTT il. Figures |2a,b) reveal that as the masses M are comparable for 
both models, the important differences come from A^cc — A^cs- Hence, we can conclude that it is 
not the actual shape of the azimuthon which determines H, but the convolution integrals A^cc 
and A'cs with the response function. In fact, A'cc and A^cs represent the overlaps of the nonlocal 
response of one dipole component with either itself or its 7r/2-rotated counterpart. Hence, the 
more of the non-rotational-symmetric shape is preserved by the nonlocal response, the larger 
is A'cc ~Ncs- Looking at both response functions considered here in Fourier domain [Fig.lJJc)] 
it is obvious that the degree of nonlocality and consequently the degree of spatial averaging is 
greater for the Gaussian response leading to smaller difference between A'cc and A'cs- 

4. Rotation Frequency - Numerical Analysis 

How good are our approximations of the rotating frequency D.1 Is the dependency in the struc- 
tural parameter a really approximatively a/ [a? + 1) for a given Aq? To answer these questions 
we compute the azimuthons U numerically and propagate them over some distance in z to ob- 
serve their rotation and estimate ii. One typical example for the Gaussian nonlocal model is 
shown in Fig. [3] The corresponding vortex soliton has a propagation constant Aq = 11.5 and 
mass M ~ 200. As its width a w 0.7 is smaller than the width of the response function we 
are definitely in the nonlocal regime. We can clearly see the reasonable agreement between 
approximate and fully numerical results [compare with Fig.fljc)]. We also tested azimuthons 
with larger masses and therefore corresponding to higher degrees of nonlocality. For instance, 
for Ao = 247, M = 2000 and a = 0.95 we find ii — 0.54 from direct numerical simulations, 
in excellent agreement with Q. = 0.52 from the approximate model [see also Fig. [He)]. Fig- 
ure|4]shows another example, this time for the Bessel nonlocal model. Again, the agreement is 
reasonable. Here, we have to use higher Ao and larger mass in order to ensure the stability of 
the vortex soliton. Like in the previous example, we operate in a nonlocal regime (a = 0.17). 
To summarize, our simple approximate model for Q., which involves nothing but a variational 
solution for the vortex soliton, allows us to predict the correct range and sign of Q. for a family 
of dipole azimuthons with given mass. Moreover, the dependency on the structural parameter 
a is also correctly modeled. 

However, there is still serious discrepancy in the values of ii|a=i. One could think of two 
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Fig. 3. Gaussian nonlocal model, M = 200, Ao = 11.5: Solid line shows the result of 
Eq. Jilt computed from the variational vortex approximation. Crosses depict results ob- 
tained from direct numerical solutions. Dashed lines represent results of the linear pertur- 
bation analysis of dipole and vortex soliton, respectively (see Sec.|5]l. 




Fig. 4. Bessel nonlocal model, M = 2000, Ao = 418: Same coding as in Fig.[3] 



possible reasons. Firstly, it could be the deviations between the actual vortex soliton and its 
variational approximation. However, it turns out that 2(A'cc ~Ncs)/M computed from the nu- 
merically obtained vortex profile is almost constant (within a few percents). Hence it has no 
effect on the overall error in D.\a=i- However, postulating the specific ansatz Eq. (|6]) we im- 
plicitly assumed a certain shape of the deformation to the soliton profile while going over from 
the vortex to azimuthons, in the limit a —> 1. Therefore it has to be this shape of vortex de- 
formation which determines the rotation frequency ii, since a vortex formally allows for all 
possible rotations (see discussion on shifting X in end of Sec. [2]i. In the next section we will 
therefore discuss the formation of azimuthons by considering it as a process of bifurcation from 
the stationary non-rotating soliton solutions. 

5. Internal Modes and Azimuthons 

Let us start with the dipole soliton D(r, 0), because here the resulting linear problem is sim- 
pler than in the vortex case. We can directly use the rotation frequency Q. as the linearization 
parameter, since for the dipole we have iio = 0. We use the following ansatz 



f/(r, ^ - Q.z) = D{r, <j)-D.z) + i^d{r, <j)-Q.z), 



(15) 



where D and d are real functions representing soliton and its small deformation, respectively. 
We then plug this ansatz into Eqs. ([T]), (|5]l and (|6]l and linearize with respect to ii. The resulting 
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Fig. 5. Dipole D and deformation d for tlie Gaussian nonlocal model, M = 200, A/j = 1 1 .8. 
The deformation d has a dipole shape, rotated by Tt/2 with respect to D. We can construct 
an approximate azimuthon U = D + i£ld for small rotation frequencies Q.. The amplitude 
ratio maxii/maxD gives the slope of Sl{a) at a = 0. 



equation of order 6'[QP) is fulfilled by the dipole soliton D. The terms of the order ff{'D} ) lead 
to the following linear equation for the deformation d: 



■Xnd- 



~A±d+ II R{\r- 



r\)D'{/)d'r'd^—D 
dip 



(16) 



In deriving the above equation we made use of the fact that the propagation constant A of the 
azimuthon depends on il as A = Ad + i^(fi^), where Ad is the propagation constant of the 
dipole soHton D. This can be seen easily from Eq. (jTa]) and M = Md + ^(ii"), L- = + ^(i2), 
/ = /d + ^(n^), and N = Nu + (?{^^), where the index "D" indicates the respective integral 
for the dipole soliton D. 

The solution ofEq.lfTSIl for the Gaussian nonlocal model and the dipole soliton with M = 200, 
and Ad = 11.8 is shown in Figure |5] The left plot shows the dipole soliton D{x,y) while the 
right shows its deformation d{x,y). As we can see from the right subplot, the deformation 
d has a shape of a dipole as well, but rotated by n/2 with respect to D. The width of d is 
about 0.85 that of D. Hence, for small rotation frequency Q. the dipole azimuthon consists 
of two orthogonally oriented dipoles, with the the relative phase shift of n/2 and amplitude 
ratio a = Q.ma\d/ma\D. In effect, the whole rotating structure is slightly elliptic, with the 
principal axis given by D. In our example, we find maxd / maxD = 1.53, hence Q. = 1.53a, 
which agrees perfectly with exact numerical results obtained from the solution of the original 
nonlinear problem Eq.([T| for a ^ (see Fig.|3]l. 

We obtain similar results for our second example, the dipole soliton in the Bessel nonlocal 
model with M = 2000, Xd = 444 (see Fig. |6]l. The ratio of the widths of a dipole D and its 
deformation d is 0.76, and Q. = 62.1a. Again, the agreement with the full numerical results of 
Fig. His excellent. 

Let us now look at the azimuthon originating (bifurcating) from the vortex soliton. For this 
purpose, we recall the Eigenvalue problem for the internal modes of the nonlinear potential 9 
which is usually treated in the context of Unear stability of nonlinear (soliton) solutions. We 
introduce a small perturbation 8V to the vortex soliton V, and plug 



\(f={V + 8V)e' 



i<j>+iAfiz 



(17) 
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Fig. 6. Same as Fig.|5]but for the Bessel nonlocal model, M = 2000, ^o = AAA. 



into Eqs. ([T]i and (|5]l and linearize with respect to the perturbation. Please note that the pertur- 
bation 5y (r, ^,z) is complex, whereas the vortex profile V{r) can be chosen real. The resulting 
evolution equation for the perturbation 5V is given by 



^^HiM^'hll"''-'^'^-''''"' 



+V / / /?(|r - /DV (r') [5V(7',z) + 8V*{7\z)\ d^ = 0. 



d , Id 

'3--A0 + -3- 
az r or 



With the ansatz 

we then derive the Eigenvalue problem for the internal modes 



5V 



dV2 



= K 



fSVi 
\8V2 



where [J 



5Vi 

5^2 



VjjR{\r-r'\)Vir')[dViir') + dV2ir')]cos[m{(l)-(t)')]d^r' 
-V JJ R{\r~r'\)V{r')[5V2{r') + 5Vi{r')]cos[m{^ - ^')]d^r' 




(18) 



(19) 



(20) 



(21) 



Real Eigenvalues of Eq. ( |20] | (ic = fc*) are termed orbitally stable and the corresponding Eigen- 
vector (SVi , 5V2) can be chosen as real function. If we perturb the vortex V with an orbitally 
stable Eigenvector, the resulting wave-function y/ can be written in the form of Eq. ^ with 
Q. = —K/m and A = Ao — K/m. Hence, we expect to find an orbitally stable internal mode with 
m = 2 corresponding to the azimuthon, where the Eigenvalue K gives n|a=i = —k/2. 

Indeed, when we solve problem ( l20l i for our two test cases numerically we find Eigenvalues 
at the expected positions, K = — 1 .4 and K = —62 respectively. The resulting estimates for 



' Please not that since | r — r' | = i/r^ + r'^ — 2rr'cos{(|) — <p'), all integrals in I2U are independent of 1 




Fig. 7. Gaussian nonlocal model, M = 200, Aq = 11.5: Solid line shows the component 
SV\, dashed line the component SV2, for ni = 2 and Eigenvalue K" = — 1.4 computed from 
Eq. J20t . The movie ( l4MBI or ll5 MBt shows the evolution of intensity and phase of the 
vortex perturbed with the Eigenvector shown above, f/|-=() = Vexp{i(j)) + SVi exp{i3(j)) + 
SV2exp{—i(j)). The chosen amplitudes of the perturbation are such that a = 0.85, the rotat- 
ing frequency is close to 0.\a^i = 0.7. 



£l\a=i fit perfectly with values obtained from direct simulations (see Figs.[3]and|4]i. The radial 
profiles of the corresponding Eigenvectors are shown in Figs. |7] and [8]). The total vorticity of 
the component ~ SVi is 2 + 1 = 3, the one of '^ dV2 is —2 +1 = —1. These total vorticities 
manifest themselves in the shapes of the radial profiles near the origin r = 0. 

As for the perturbation of the stationary dipole, it is possible to construct azimuthons in the 
vicinity of the vortex (awl) from dV: 



{/(r, (/))|,=o = [V{r) + dVi (r)e""^ + 5y2(r)e-""^] e 



im^l i^ 



(22) 



Used as an initial condition to the propagation equation this object is expected to rotate with 
£l\a=i = —K/m. Unlike in the previous case of the stationary dipole, the amplitude of the per- 
turbation is not fixecQ but will eventually determine the value of the structural parameter a. 
Generally speaking, the smaller the resulting a the greater the error in the constructed ini- 
tial condition. However, the great robustness of the azimuthons allows one to use the initial 
condition ( l22b for quite large perturbation amplitudes (resulting a ^ 0.5). Those "bad" initial 



-Just the ratio between the components SV[ and SVi is fixed. 
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Fig. 8. Bessel nonlocal model, M = 2000, Aq = 418: Same as in Fig.|7l m = 2,K = -62. 
The movie ( 4MB or 15 MB ) shows the evolution of the perturbed vortex with a = 0.9, the 
rotating frequency is close to ri|a=i =31. 



conditions result in oscillations of the azimuthon upon propagation. However, the azimuthon 
is structurally stable and does not decay into other solutions like the single-hump ground state. 
Moreover, such initial conditions play a role of excellent "initial guesses" for solver routines to 
find numerically exact azimuthons. In the movies shown in Figs. [T] and |8] we choose moderate 
amplitudes of the perturbations in order to get almost perfect azimuthons. The corresponding 
structural parameters a > 0.85 also guarantee that the observable rotating frequencies are very 
close to their limiting values Q.\a=i- 

The results of the linear stability analysis of the vortex soliton offer an easy explanation for 
the observed systematic deviation of the estimation for n|£,=i by Eq. ( fTTT i from its true value. 
As mentioned before, the shape of the perturbation to the vortex soliton determines n|a=i. In 
our ansatz (|9]l we only account for a perturbation similar to the component '^ 5V2 and "forget" 
the other component ~ SVi with total vorticity 3, which leads to the observed deviation. The 
larger the "forgotten" component, the larger is this deviation (compare Figs. |7] [3] and Figs.[8l 

a. 

There is more to say about the internal modes of the vortex soliton. On one hand, solving the 
Eigenvalue problem ( l20l i for other vorticities m offers an exhaustive method for finding families 
of azimuthons which originate from a vortex soliton. For example, with our method it was 
straightforward to identify a rotating triple-hump azimuthon (m = 3, K" = 3.8), which rotates 
with Q. = —K/m = —1.27. This solution bifurcate from the M = 200 vortex in the Gaussian 




Fig. 9. Gaussian nonlocal model, M = 200, Aq = 11. 5: Same as in Fig.[7] m = 3, K = 3.8. 
The movie (4 MB or 15 MB I shows the evolution of perturbed vortex C/|;=o = V exp{i(j>) + 
SVi exp((4(|)) + 5V2exp(— (20), a = 0.9, the rotating frequency isclose tofi|a=i = — 1.27. 



nonlocal model. The components of the generating Eigenvector are shown in Fig. |9l Again, we 
can easily construct various azimuthons with a « 1 from this Eigenvector The movie in Fig.|9] 
shows an example of the propagation of the triple-hump azimuthon with a ~ 0.9. 

It should be stressed that not all orbitally stable Eigenvalues can be linked to a family of 
azimuthons. This is obvious for Eigenvalues with jfcj > |Ao| in the continuous part of the spec- 
trum. Hence, we can conclude that |ii|a=i | < |Ao|/m for azimuthons bifurcating from a single- 
charged vortejo. However, bound orbitally stable Eigenvalues do not necessary indicate an 
emanating family of azimuthons. For example, the M = 200 vortex in the Gaussian nonlocal 
model possesses several other bounded internal modes for m = 2, e.g. JC = 7. If we perturb 
the vortex with the corresponding Eigenvector the resulting structure is a double-hump rotating 
with n = — 3.5. Closer inspection reveals that this structure decays and hence does not belong to 
another family of dipole azimuthons. Nevertheless, since those decaying rotating structures can 
survive over large propagation distances they may be indistinguishable from true azimuthons 
in the experimental conditions. This issue is currently under detailed investigation. 



Please note that the parameter m determines the number of humps of the rotating structure. 



6. Conclusion 

In conclusion, we analyzed properties of azimuthons, i.e. localized rotating nonlinear waves in 
nonlocal nonlinear media. We showed that the frequency of angular rotation crucially depends 
on the nonlocal response function. Starting with the single charge vortex soliton and by em- 
ploying a simple variational ansatz we were able to predict accurately the frequencies of the 
rotating dipole azimuthon. In addition, we could explain how the actual shape of the response 
function determines the rotation frequencies. Further on, we computed exact dipole azimuthon 
solutions and their rotation frequencies numerically and showed that in the limits of maximal 
and/or minimal azimuthal amplitude modulation, i.e., close to the dipole or vortex soliton, the 
rotation frequency is determined uniquely by Eigenvalues of the bound modes of the linearized 
version of the respective stationary nonlocal solution. Moreover, the intensity profile of the 
resulting azimuthon can be constructed from the corresponding linear Eigen-solution. This of- 
fers a straightforward and exhaustive method to identify rotating soliton solutions in a given 
nonlinear medium. 
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